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We present a study of the influence of different types of disorder on systems in the Ising uni- 
versality class by employing both a dynamical field theory approach and extensive Monte Carlo 
simulations. We reproduce some well known results for the case of quenched disorder (random 
temperature and random field), and analyze the effect of four different types of time-dependent 
disorder scarcely studied so far in the literature. Some of them are of obvious experimental and 
theoretical relevance (as for example, globally fluctuating temperatures or random fields). All the 
predictions coming from our field theoretical analysis are fully confirmed by extensive simulations 
in two and three dimensions, and novel qualitatively different, non-Ising transitions are reported. 
Possible experimental setups designed to explore the described phenomenologies are also briefly 
discussed. 

PACS numbers: 05.10 Cc, 05.10 Ln, 05.70 Ph, 64.60 Cn, 75.10 HK 



The Ising model has played a central role in the devel- 
opment of modern statistical mechanics. Its universality 
class includes a broad variety of models and real systems, 
as magnetic and reaction-diffusion systems, spin glasses, 
neural networks, coupled maps, and protein folding mod- 
els jj],U . The analysis of the robustness of this universal- 
ity class upon the introduction of disorder is of outmost 
importance in many different contexts. In fact, there 
are plenty of studies devoted to determining whether, 
and in which way, quenched disorder is a relevant per- 
turbation at the Ising renormalization group (RG) fixed 
point While a lot of attention has been paid to 

the study of quenched-in-space disorder Q|| in its dif- 
ferent variants, quenched-in-time disorder (i.e. homoge- 
neous disorder changing randomly in time) has not yet 
been thoroughly studied (some exceptions to this asser- 
tion are i see a l so || ) ■ This is surprising as stochasti- 
cally time-changing temperatures, or magnetic fields, are 
quite straightforward idealizations of realistic situations. 

Here we study the effect of different types of disorder 
on the Ising universality class by using field theoretical 
tools and extensive Monte Carlo simulations. Quenched- 
in-space disorder will be considered for the sake of com- 
pleteness (in particular, quenched-in-space temperature 
(QST), and quenched-in-space random field (QSF)), but 
we will chiefly focus on the novel and rich phenomenol- 
ogy of temporally disordered Ising models, and will 
analyze extensively the following four types of disor- 
der: (1) Qucnchcd-in-timc, uniform-in-space, random 
temperature (QTT). (2) Quenched- in-space-time random 
temperature (QSTT). (3) Quenched-in-time, uniform-in- 
space, random (magnetic) field (QTF). (4) Quenched-in- 
space-time random (magnetic) field (QSTF). Physically, 
the difference between the odd and the even cases stems 
from the fact that in the first ones the fluctuating exter- 
nal parameter couples globally to the system, while even 
cases represent situations in which the fluctuating pa- 



rameter is inhomogeneous (changing in space and time). 
In terms on Monte Carlo implementations, this difference 
can be enviewed as associated with two different types of 
updating: if a value of the random variable is selected for 
each spin update (sequential updating), one is placed nat- 
urally at the second (or fourth) case, while if the random 
variable is extracted after updating of the whole lattice 
(this is the case, for example for parallel updating) the 
situation is described by the first (or third). In this way 
the type of updating can indirectly determine the nature 
of the disorder, and eventually the emerging behavior. 

Our main conclusions are (listed in table 1): (i) QTT 
is a relevant perturbation for d > 2, but irrelevant in 
d = 2. We find strong numerical evidence of continu- 
ously varying exponents as a function of the width of the 
temperature distribution. A qualitatively different phe- 
nomenology is found for strong disorder (i.e. for widths 
above a certain threshold), (ii) QSTT is always an irrel- 
evant perturbation, (iii) QTF is relevant for d > 1; and 
a phenomenology completely different from the standard 
Ising one emerges, (iv) QSTF is marginal for d = 4, and 
irrelevant otherwise. 

Let us start discussing the disordered-temperature 
cases. The most basic result for QST disorder is the 
Harris criterion || that can be sketched as follows: If T 
is a random uncorrelated || quenched-in-space variable, 
the amplitude of its fluctuations per unit volume, /, in 
a region of linear size £ is proportional to £ _d / 2 . At a 
distance A of the critical point and taking £ to be the 
correlation length, one has £ - and / oc A~ . If 

vd < 2 temperature fluctuations control the critical scal- 
ing, while they are irrelevant otherwise fiofl . To study the 
QTT disorder, one may replace £ by a characteristic time 
to ~ A 1/2 and follow analogous steps to obtain vz < 2 as 
the proper relevance criterion. 

In order to make more rigorous and general the previ- 
ous arguments let us now re-derive them in a field the- 
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oretical context M. First we cast the Ising model into 
a Ginzburg-Landau- Wilson (GLW) dynamical functional 
[fill, whose associated action is 



over a magnetic-field Gaussian distribution (centered at 
zero and of width 5) ||, we obtain the terms to be 
added to Eq. (g): (i) 5 f <bc[J dtijj(x,t)] 2 for QSF (ii) 



S = J dxdt{^ip 2 (x,t) - ip(x, t)[0(x, t)-T(f>- 
^V 2 0(x,t) + |0 3 (x,i)+^(x,i)]} 



where <$> is the magnetization field, ip is the response func- 
tion, A and g are parameters, r oc T — T c includes the 
relevant dependence on the temperature |T^| , and h is a 
magnetic field. 

In order to include the effect of a disordered tempera- 
ture we substitute r by a function t(x, t). Then we can 
resort to one amongtwo alternatives: one is to employ 
the "replica" trick Q, while here we follow the equiv- 
alent dynamical approach introduced some time ago by 
De Dominicis fl3]| . In fact, the dynamical approach is the 
most natural one for the study of time-dependent disor- 
der. The next step is to average the dynamical generating 
functional associated with Eq. (g) over the distribution 
of disorder. In particular, writing r(x, t) — tq + 6t(x, t) 
(where To is the averaged temperature and 5r(x, t) is 
a zero-mean, a variance, Gaussian distributed random 
variable and performing the Gaussian integration, 
one obtains an effective action identical to Eq. (g) except 
for the replacement r — > To and for the presence of the fol- 
lowing additional disorder-dependent (nonlocal) terms: 



g J dx 

a dt 



dt(j)(x, t)ip(x, t) 
dx<p(x, t)ip{x, t) 



for QST (2) 



for QTT (3) 



dt[<P(x,t)i;(x,t)} 2 for QSTT. 



(4) 



Let us analyze separately these three cases. For QST 
the naive dimension of the coefficient of the extra (non- 
Markovian) term is such that Eq. (J2j> is dimension- 
less. As J dt J dXTQ(fnl> has also to be dimensionless, 
[/ dt(j)tp} = A 4 - 1 /". Substituting this into the previous 
expression one obtains [a] = A 2 /"~ d , and therefore a is 
relevant (in the infrared limit ]T^]) if dv < 2 (reproducing 
the Harris criterion). Proceeding analogously for QTT 
one obtains vz < 2 and for QSTT v{d + z) < 2 as rel- 
evance criteria. Substituting the known Ising exponent 
values |l2], one finds that: (see table) QST is a marginal 
perturbation at d = 2 jl4j and relevant at d — 3, as al- 
ready known QTT is relevant for d > 3, while it 
is irrelevant in d = 2. QSTT is irrelevant in any dimen- 
sion. This prediction implies that homogeneous globally 
changing temperatures (QTT) are far more effective in 
altering Ising behavior than locally changing tempera- 
tures (QSTT). 

For the case of random fields, we proceed along anal- 
ogous lines. By averaging the generating functional 



Sjdt[J dx^(x, t)] for QTF, and (iii) SjdxJ dt^ 2 (x, t) 
and for QSTF. 

By noticing that the dimension of J dtip(pc,t) at the 
(1) pure Ising RG fixed point in terms of A is | — 3- (where 
7 is the susceptibility exponent), it is easy to derive the 
following relevance criteria: j/v > for QSF, dv+~/ > vz 
for QTF, and 7 > vz for QSTF. Substituting the Ising 
exponents [|2| we are lead to the following predictions 
(see table 1): (i) QSF is relevant for d > 1. For this case, 
the Imry-Ma argument establishes that d = 2 is the lower 
critical dimension of the QSF universality class fll|). (ii) 
QTF is relevant in all dimensions. Some previous results 
about this model can be found [^7); in particular, it is 
expected to exhibit a transition even in d — 1 , and a non 
Ising-like phase transition has been reported in d = 2 
0. (iii) QSTF is marginal only at d = 4, and irrelevant 
otherwise (in fact, it just changes the noise amplitude, 
therefore it does not introduce any new operator). Anal- 
ogously to what happens in the case of disordered tem- 
peratures, homogeneous global magnetic fields (QTF) are 
more effective in affecting the universal properties of the 
pure Ising class than magnetic fields varying both in time 
and in space (QSTF). 

Let us finally stress that all the relevance criteria are 
determined using anomalous exponents, therefore at the 
Ising and not at the Gaussian fixed point. By perform- 
ing naive power counting one can easily see, that QTT 
and QTF are also relevant perturbations at the Gaus- 
sian, mean-field, fixed point. Hence they are relevant 
even above the Ising upper critical dimension, d c = 4. 

Having established all these relevance criteria, the next 
natural step would be to determine the new upper criti- 
cal dimensions, d Cl in the cases where disorder is relevant, 
and perform an e— expansion in each case around the cor- 
responding d c . For example, for QSF the extra nonlocal 
operator scales as A 2 at d — 4 (it is more relevant than 
the Ising nonlincarities), and accordingly the upper crit- 
ical dimension is shifted from 4 to 6. To every order in 
an e— expansion the exponents of the QSF model in d 
dimensions are identical to those of the pure problem in 
d — 2 dimensions 17 J] . 

Our attempts to perform similar detailed RG calcula- 
tions in the other cases have turned out to be extremely 
hard and, so far, inconclusive [Q. For instance, for QTT 
and QTF the degree of divergence of the nonlocal opera- 
tor is insensitive to dimensionality above the Ising critical 
dimension d c — 4. Consequently there is no upper critical 
dimension for these types of disorder which are strongly 
relevant even at high dimensions. Hence, even finding a 
sound mean-field approach, around which one could per- 
form some sort of perturbative approach, turns out to be 
a nontrivial task. The only remaining possibility we can 
think of would be to perform a double expansion anal- 
ogous to that in In any case, at the present stage 
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of our research, we can just predict whether disorder is 
relevant or not, but cannot make sensible predictions con- 
cerning eventual new RG fixed points and the nature of 
new transitions. 

In order to test the above predictions, and explore 
the eventual existence of new universality classes and/or 
novel phcnomenologies, we have performed extensive 
Monte Carlo simulations combined with finite size scal- 
ing analysis in d — 2 and d = 3 with different types of 
disorder. We use: system sizes up to L = 512 (L = 90) in 
d = 2 (d = 3), periodic boundary conditions, Metropolis 
H transition rates, and run up to 3 x 10 7 Monte Carlo 
steps. Here we report the main results, further details 
and discussions will be presented elsewhere Q: 

(i) For disordered temperature cases we have per- 
formed simulations using both QTT and QSTT. The 
temperature is distributed homogeneously and symmet- 
rically around a central (tunable) value To and width <5 
||. For both types of disorder in d = 2 we obtain ex- 
ponents compatible with the respective systems being in 
the Ising class; see inset in Figure 1, where we study 
the scaling of the susceptibility peak (maximum value of 
the magnetization fluctuations) with system size L. The 
model studied in |H| can also be classified as QSTT in 
d = 2, and in fact, it was reported to exhibit Ising scal- 
ing. In d = 3 for QSTT from the curve ^vsl we observe 
very good scaling using the three-dimensional ^jv Ising 
value. On the other hand, for QTT we observe exponents 
varying continously with 5, changing from the Ising value 
"f/v = 1.97(1) for 5 = 0, to j/v = 3 for values of 8 « 1.5. 
There is no evidence of crossovers and the variation in ex- 
ponent values is rather broad. Therefore we believe that 
these are the real asymptotic values and that accord- 
ingly there should be a line of RG fixed points. Caution 
should be used before concluding that the exponents ac- 
tually change continuously. Corrections to scaling could 
be responsible for such a variation j22|. For values of <5 
larger than a certain threshold around 1.5 the exponent 
7/1/ saturates to 3 implying that the magnetization his- 
tograms do not converge to delta Dirac functions in the 
thermodynamic limit: they remain broad, corresponding 
to external-fluctuations-controlled transitions. 

(ii) For the case of random fields, QTF is expected 
to be non-Ising in d = 2, while QSTF disorder is pre- 
dicted to be an irrelevant perturbation. In fact, for QSTF 
our simulations show a perfect agreement with Ising ex- 
ponents for all the measured magnitudes; non-universal 
properties depend on the disorder width, 5, but criti- 
cal exponents remain unchanged. The models studied 
in [p3|2]| can also be ascribed to this class of disorder, 
and were found to be Ising like in d = 2 in accordance 
with our prediction. On the other hand for QTF we 
observe a completely different phenomenology: There is 
a continuous phase transition whose features are unam- 
biguously different from the corresponding Ising ones. In 
particular, it is a sort of dynamical transition controlled 
by the fluctuations of the randomly changing magnetic 
field. For example, in figure 2, we show the dependence 



of the susceptibility per unit surface as a function of sys- 
tem size, for different values of S. After a finite-size-effect 
transient, the curves converge to a fixed value, indicating 
the existence of broad magnetization histograms in the 
thermodynamic limit, and implying j/v = 2. Still, mag- 
nitudes as the specific heat or the susceptibility show sin- 
gular peaks at the critical point, and there is a dynamical 
spontaneous symmetry breaking below the critical tem- 
perature 20 24]]. We have reproduced qualitatively this 
QTF phenomenology by means of a dynamical mean-field 
approach to be reported elsewhere [^0| , and observed that 
the nature of this transition does not seem to be very sen- 
sible to dimensionality changes. 

Summing up: we have presented a rather simple gen- 
eral field theoretical scheme allowing us to predict the rel- 
evance of different types of disorder in the Ising universal- 
ity class. Some of the conclusions where already known, 
and provide a validity test for our approach. Using this 
scheme we have studied the effect of different types of 
time-dependent disorder in the Ising class. These disor- 
dered models describe rather natural situations in which 
systems are subjected to the influence of some globally 
(or locally) fluctuating external agent (randomly vary- 
ing temperatures or magnetic fields). All our predictions 
are fully confirmed by means of extensive Monte Carlo 
simulations; and novel quite interesting phenomenologies 
appear for the case of homogeneously fluctuating tem- 
peratures and magnetic fields. 

In order to experimentally verify our predictions one 
could take, for example, vacuum deposited ultrathin Fe 
films on W(110) (i.e. a nearly ideal two-dimensional mag- 
netic Ising system p5[|), and in the same way that this 
material has been set to study hysteresis in the presence 
of an oscillating magnetic field pq] , it can be analyzed 
under the effect of a randomly varying magnetic field 
(generated, for example, by means of an electric circuit 
in a chaotic regime) . A detailed comparison between the 
characteristic relaxation time-scale, and the typical mag- 
netic field variation scale is of outmost importance in 
order to determine the applicability of our results to this 
type of experiments. 
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log L 

FIG. 1. Susceptibility at the critical point as a function 
of system size for disordered-temperature models in d = 3 
(d = 2 in the inset). QSTT exhibits Ising scaling in d = 3. 
For QTT the slopes (fixing ~y/v) change continuously with 
disorder amplitude in d — 3, but not in d = 2 (see inset). 
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FIG. 2. Susceptibility per unit surface as a function of L 
(in log-log scale) for QTF disorder. For large sizes the curves 
converge to a constant indicating the existence of broad mag- 
netization histograms in the thermodynamic limit. 
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TABLE I. Criteria for the relevance of different types of 
disorder in the Ising universality class. The first row in each 
case specifies the general criteria (equalities correspond to 
marginal cases), while in the second we specify the dimen- 
sions at which such criteria are fulfilled. 
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